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Abstract 

Arbitrary amplitude nonlinear ion waves is investigated in an extremely degenerate electron- 
positron- ion plasma with relativistic electrons/positrons and dynamic cold ions using Sagdeev 
pseudo-potential method in framework of quantum hydrodynamics model. The matching criteria 
of existence for ion solitary as well as periodic nonlinear excitations is studied numerically in 
terms of the relativistic degeneracy parameter, relative positron concentration and the relativistic 
normalized propagation-speed range of these waves are obtained. It is remarked that the electron 
relativistic degeneracy and relative positron concentration in such plasmas have significant effects 
on nonlinear wave propagations in relativistically degenerate Fermi-Dirac plasma. Our results are 
in good agreement with the previously reported ones obtained using semiclassical Thomas-Fermi 
approximation in dense plasmas with non-relativistic and ultra-relativistic electrons and positrons. 
Current findings can be appropriate for the study of astrophysical superdense compact objects such 
as white dwarfs. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. INTRODUCTION 



Ion waves are one of the most prominent features of ionized environments with abun- 
dant occurrence in astrophysical entities and laboratory plasmas. The characterization 
of such waves began early in 1961 by work of Vedenov et.al PJ using pseudo-potential 
method. Another approach which is widely employed to investigate the asymptotic behav- 
ior of nonlinear excitations is the so-called reductive perturbation technique 2|-|6] . However, 
the later approach is suitable only for the study of small-amplitude nonlinear perturbations 
and is usually applied for treatment of plasma waves in a state very close to thermody- 
namic equilibrium. In order to achieve satisfactory results applicable to real environments 
using perturbation method, it is needed to consider the higher-order approximations for 
dependent plasma variables |z-9]. In recent years there has been many investigations on 
arbitrary-amplitude nonlinear ion solitary (ISWs), ion periodic (IPWs) as well as electro- 
static solitary waves (ESWs) in different plasma environments 1CH14I]. Furthermore 
small amplitude propagation as well as interaction of ISWs has been considered in Ref. [1| 
in magnetized pair plasmas with the electrons and positrons possessing either non-relativistic 
or ultra-relativistic degeneracy pressure. It is however the aim of current work to study the 
propagation of ion waves in dense plasmas with arbitrary degree of relativistic degeneracy 
in three component electron-positron-ion dense plasmas. 

Because of accommodation of variety of nonlinear excitations in pair plasmas with com- 
ponents of equal mass and charge, and also because of their inevitable applications in astro- 
physics, the electron-positron- ion plasma has attracted outstanding attention in recent years 
16N20|. For instance, electron-positron-ion plasmas exist in active galactic nuclei 21], pul- 
sar magnetospheres 22] and many dense astronomical objects like neutron stars and white 



dwarfs 
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24] . It is also believed that the electron-positron-ion plasma may play an impor 



tant role in the beginning and evolution of Universe 
and studied in laboratory 
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251 ] . Pair plasmas have been produced 



fined in a Tokamak 



30, 



-|29|. The electron-positron-ion plasma can also be easily con- 
31] . Calculations 32] reveal that even in dense electron-positron- ion 
plasmas ion waves can have significant lifetime to propagate in such environments, despite 
the high rate of annihilation processes. 



On the other hand, in compact dense astrophysical objects 



33] and in ordinary metals 



and semiconductors [34| , fermion constituents of plasma are under extreme condition called 
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the degeneracy state. In a degenerated plasma the quantum and relativistic effects become 
important and the much higher degenera cy p ressure which arises due to the Pauli-exclusion 



mechanism should be taken into account 



351 ] . The quantum effects become more important, 



specifically, when the inter-particle distances become less than the characteristic de Broglie 
thermal wavelength, = h/ftnmkBT) 1 / 2 36|. Chandrasekhar [37|, using the Fermi-Dirac 
statistics for electrons, has given a mathematical criteria under which a white dwarf can 
be considered as completely degenerate. In the zero-temperature Fermi-gas, used for a 
completely degenerate plasma, it is assumed that the thermal temperature of plasma species 
are , although finite, but fall much below the well defined characteristic Fermi-temperature 
. Although, the zero-temperature model may seem a simplifying assumption, however, it 



hold for many dense astrophysical entities such as white dwarfs [39|, |40| . It should be noted 
that to have a more realistic model accounting for the electron/positron temperature, one 
has to revisit the finite temperature Fermi-Dirac model which can be our future work. 

It has been shown that, the electron degeneracy in massive white dwarfs which opposes the 
catastrophic gravitational inward pressure may give rise to softened relativistic degeneracy 
pressure leading to the ultimate collapse of massive star in some conditions [39]. There are 
some studies, using quantum hydrodynamics, which deal with astrophysical models consid- 



ering only the non-relativistic electron degeneracy case |4lH44|. Whereas, in a typical white 
dwarf with the electron density as high as 10 28 cm~ 3 (p ~ 10 7 (gr/cm 3 )), the degeneracy 
pressure is extreme and both quantum and relativistic considerations have to be invoked. 
Application of semiclassical Thomas-Fermi approximation indicates that the dynamics of 
nonlinear ion excitations in the electron-positron-ion plasma can be quite different when the 
electrons/positrons are under ultra- relativistic degeneracy pressure from that under normal 



dynamics model 
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degeneracy pressure 45]. Previously reported investigations, based on the quantum hydro- 



49j, consider normal degeneracy case encountered in ordinary metals 



and semiconductors. In this study, we investigate a degenerated plasma for a whole range of 
the relativity parameter (degree of electron degeneracy) which can extend the applications 
of quantum hydrodynamics model used in mentioned references to extreme conditions found 
in compact astrophysical objects. We also use the Fermi-Dirac statistics for electrons and 
positrons instead of the simplified Thomas-Fermi approximation. The presentation of the 
article is as follows. The basic normalized plasma equations are introduced in section [III 
Nonlinear arbitrary-amplitude solutions are derived in section IIHI Numerical analysis and 
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discussions is given in section [IV] and final remarks are presented in section |VJ 



II. QUANTUM HYDRODYNAMICS MODEL 

In this section we use quantum hydrodynamics (QHD) model to study the nonlinear 
dynamics in a superdense completely degenerate electron-positron-ion plasma. The elec- 
trons/positrons are considered relativistic and follow the Fermi-Dirac statistics and the 
positively charged cold-ions are dynamic. In this model the rate of electron(positron)-ion 
collisions are limited due to Pauli-blocking process and the pair-annihilations can be ignored 
32I . We therefore write the conventional QHD fluid equations in dimensional form as 
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where, n it j, and v it j denote the ions and j-species number- density, mass and velocities 
with % for ions and j = {e,p} with Sj = {— 1, +1}, for electrons and positrons, respectively. 
Also 7, Tj and H, respectively, present the adiabatic constant, ion-temperature and the scaled 
Plank-constant. The normalized set of QHD equations follow from the following scaling 

x y x, t y , Tt e pi y Ti e p jji e Q 1 Vi^ e ^p y Vi^ e ^pC STl <p y <p . (2j 

where, c, u p i = -\/47re 2 n e o/mj and c sr = \J m e c 2 /rrii are the vacuum light speed and the 
characteristic plasma frequency and sound-speed, respectively, and the parameter n e0 de- 
notes the electron equilibrium density. The bar symbols denote the dimensionless quantities 
which are omitted in forthcoming algebra for simplicity. In a relativistic completely degener- 
ate Fermi-Dirac plasma, the electron/positron relativistic degeneracy pressure is expressed 



in the following form 



37] 
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TXVTQC 



\Rj [2R] - 3) yjl + R) + Ssinh" 1 ^ } , (3) 

in which Rj = (nj/no) 1 ^ 3 (no = 8?r ^g c ~ 5.9 x 10 29 cm~ 3 ) is the ratio of electron/positron 
Fermi relativistic momentum to m e c. In the limits of very small and very large values of Rj 



we obtain 



2 



Pi={ *Wki W^°)>. (4 ) 

Thus, the normalized QHD set of plasma equations including the electron/poitron tunneling 



effect by making use of the fact that (l/rij)d x Pj = d x yjl + i? 2 , Rj = R^ajUj) 1 / 3 and 
Oij = {l,ot} for equilibrium number- densities of electrons and positrons, respectively, with 
a — n p0 /n e0 being the relative positron to electron number-density, can be written in the 
following form 
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where, we have introduced the normalized degeneracy parameter, R = (neo/no) 1 / 3 , as a 
measure of the relativistic degeneracy effects i.e. so called the normalized relativity parame- 
ter. Therefore, the parameter R^ denotes the relative electron number density with respect 
to the base value n . 

Furthermore, the parameter, (m e /mi)H^ /2 = (m e /mj)/i 2 a;pj/(m e c 2 ) 2 , is the quantum 
diffraction coefficient which is related to the relativity parameter via; {m e /m,i)H^ jl ~ 4.6 x 
10~ 10 i?Q. It should be noted that the parameter Rq is related also to the mass-density (of 
white dwarf, for instance) through the electron number-density with relation p ~ 2m p n e o(l — 
a) or p(gr/cm 3 ) = (1 — a)po-R[j with Po{gr/cm 3 ) ~ 1.97 x 10 6 , where, m v is the proton mass. 
The density p is i n the range of a mass-density of a typical white dwarf and a contribution 
of (1 — a) has been included in mass-density definition to account for the electron-positron 
pair production/anihilation. The density of typical white dwarfs can be in the range 10 5 < 
p(gr/cm 3 ) < 10 8 , which, neglecting electron-positron pair production/anihilation (a <^ 1), 
results in values of 0.37 < Ro < 8 for the relativity parameter of a typical white dwarf. S. 
Chandrasekhar, in ground of Fermi-Dirac statistics, has proven that for a white dwarf with a 
mass density p, the electron degeneracy pressure turns from P e oc p 5 ^ 3 dependence for normal 
degeneracy (R < 1) to P e oc p 4 / 3 dependence for ultra-relativistic electron degeneracy 
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(Ro 3> 1) 39[. In this calculation we show that this behavior has also a profound effect on 
the propagation and stability of nonlinear ion-acoustic excitations in a white dwarf. Now, 
neglecting the terms containing the small mass-ratio, m e /mj and assuming that k^Ti <C 
m e c 2 , we arrive at the following simplified set of basic normalized equations 
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In order to get the steady state solution to the hydrodynamics normalized set of 
fluid equations, Eqs. ([6]), we change to a stationary frame of soliton by changing 
the coordinate to £ = x — Mt, with, M being the constant normalized speed of 
wave train relative to the quantit y, c., r (not be confused with the Much-number 
parameter as pointed out in Ref. [501]). Therefore, after trivial integration with 



the appropriate boundary requirements; lim 



{1 — a, 1, a} and lim = 0, 

fi.e.D-»0 



we obtain 



M[n t -( 
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Solving Eqs. ((Tj) for the ion, electron and positron number-densities in terms of electrostatic 
potential leads to the following energy-density relations 
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These equations are completely different 



classical Thomas- Fermi model 



15, 



32 



44 



rom the energy relations frequently used in semi- 



45], as it should be. Note that, unlike for Thomas- 



Fermi model for plasma, in Eqs. (JSJ) there is no explicit contribution for electron/positron 
(Fermi) temperatures. In fact there is a relation between the Fermi-temperatures and the 
number-densities of quantum species i.e. Tpj = Epj/ks = (ft 2 / (2mjkB))(3n 2 nj j0 ) 3 / 2 . In the 
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following section we will find the appropriate pseudo-potential, based on relations Eqs. (JHJ) 
and Poisson's condition, to evaluate the propagation of nonlinear ion-acoustic waves in the 
plasma described by Eq. (EJ). 



III. PSEUDO-POTENTIAL APPROACH 



Proceeding with the pseudo-potential approach, we aim at obtaining the appropriate 
energy integral based on the zero-temperature Fermi-Dirac statistics for electrons and 
positrons. The solution to this integral yields a description for the propagation of arbitrary- 
amplitude ISWs as well as IPWs periodic structures in a relativistically degenerate electron- 
positron-ion plasma. Using the density relations given in Eq. flB]) and substituting into 
Poisson's relation, we obtain 

3/2 
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Now, multiplying Eq. ((9]) by and integrating once leads to the desired energy integral in 
terms of pseudo-coordinates (0, £) 
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or, the well-known energy integral reads as 



(ii) 



7 



with the long expression for the pseudo-potential U(<f>) given as 
U{<j>) = | (sM(a - l)^M 2 -20 + 3i?o 3 In M (M 2 - 20) - SR^a" 1 In M (M 2 - 20) 
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The pseudo-potential given U(4>) and its first derivative d^U{4>) both vanish at = 0, 
as expected. From density relations given in Eqs. dHJ, we may derive the possible maxi- 
mum/minimum values of the electrostatic potential in which the pseudo-particle of as unit 
mass is confined, these maximum/minimum values are 



hp 
2 



J ml — 
b m2 = ±1 + y/T 
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The possibility of solitary waves then relies on the following boundary conditions to met, 
simultaneously 
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4>=o 



For the existence of the ISWs, it is further required that for at least one either maximum 
or minimum nonzero 0-value, we have U(4> m ) = 0, so that for every value of in the range 
m > > or m < < 0, U(<p) is negative. In such a case, there will be a potential 
minimum, which describes the propagation of an ion-acoustic solitary excitation. 
Also, for IPWs excitations to propagate it is required that 



dU{(j>) 



0=0 



0. 



d 2 U((f)) 



> 0, 
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On the other hand, it is also required that for at least one either maximum or minimum 
nonzero 0-value, we have U(<f> m ) = 0, so that for every value of <ft m the range (f) m > <p > 
or (f) m < < 0, U(4>) is positive. In such a condition we will have potential maximum which 
describes the propagation of an ion-acoustic periodic structure. 

It should be realized already that the derivation of an analytic expression for the matching 
critical normalized propagation-speed, M, form the second derivative of the pseudo-potential 
given by Eq. [12] is not possible, hence, we will indulge for the numerical analysis which is 
presented in the next section. However, some information may be obtained via the small- 
amplitude limit. Expanding the potential U(<p) near cf) = 0, we obtain 

J I" TV" 

W) = + \4>\ (16) 



where, Ug = U"(<p = 0) and U(? = U'"{<j> = 0). Inserting into Eq. 
obtain (provided that Uq < 0) a solitary solution of the form (see [511 ]) 



and integrating, we 



* K » = - 3 f (17) 

This pulse profile is equivalent to the stationary soliton solution of the Korteweg-de Vries 
(KdV) equation, obtained by use of the reductive perturbation method [521 ] . It is important 



to notice that the soliton width A = 2/^/—Uq and amplitude O — ~3Uq/Uq satisfy 
4>qA 2 = 12/Uq = est., as expected from the KdV theory. 



IV. DISCUSSION 



Figure 1 shows the regions of existence of ISWs/IPWs in yellow/blue (right /left region) 
for different values of relative positron concentration, a, in M-Rq plane. It is apparent that 
the allowed matching normalized propagation-speed range, M, is significantly affected by 
the relativity parameter, Rq as well as the fractional positron concentration, a. It is observed 
that the IPWs posses lower matching normalized propagation-speed value, M, compared to 
ISWs for all values of Rq and the relative positron density, a. It is also remarked that, as 
the relativity parameter increases, the matching normalized propagation-speed range, M, for 
both nonlinear excitation types increase and the allowed matching normalized propagation- 
speed range, M, tends to higher values, while, the effect of the increase in relative positron 
density on the IPWs is to shift the matching normalized propagation-speed range, M, to 
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lower values for all relativity parameter values, while this change on the ISWs is completely 



different. The similar plots have been presented in Ref. [53] for the relativistically degen- 
erate electron-ion plasma. The increase in the value of small-a is to enhance the matching 
normalized propagation-speed values, while it is to decrease the values of matching nor- 
malized propagation-speed for large-a values. The plots shown in Fig. 2 confirm that 
there is a maximum matching normalized propagation-speed value for ISWs in M-a plane 
shown for two extreme low- and high-relativistic degeneracy values, namely, Rq = 1 and 
-Ro — 50. This maximum matching normalized propagation-speed value is M m ~ 3.75 for 
Ro = 1 while it is diminished to M m ~ 2 for Rq = 50 which correspond to ultra-relativistic 
degeneracy case. Note that, the maximum values found here roughly agree with the corre- 
sponding values reported for Thomas-Fermi electron-positron-ion plasmas in non-relativistic 
and ultra-relativistic degeneracy limits (see Fig. 1(a) in Ref. [32]). 

It is further noticed from Fig. 2 that the lower part of plots (smaller a-value) is affected 
more by change in the relativistic factor Rq. It should be mentioned that, as the relativity 
parameter, Rq, increases the quantum effects (via the parameter H r ) become dominant, 
hence, the quantum force term in (Eq. (JHJ)) can never be ignored for very large value of Rq. 
For instance, for the very large value of Rq = 10 3 (which may be realized for a collapsing 
white dwarfs or neutron stars) we get a contribution of order (m e /mi)H^/2 ~ 0.46 for the 
quantum correction coefficient. Therefore, while Fig. 2(a) qualitatively explains the effect of 
relativistic degeneracy effects on propagation of nonlinear waves in Fermi-Dirac pair-plasma, 
a correction must be made for the very highe values of Rq. However, it has been shown that 
the quantum diffraction parameter does not modify the wave amplitude (which increases 
with increase in the value of relative positron concentration, a) and makes only corrections 



to the width of the nonlinear wave 



32] . This leaves only a minor correction for quantum 



correction which effects all values of positron concentration equally, and therefore does not 
alter our current qualitative approach. 

The profiles of pseudo-potentials and their variation with respect to the 
relativity parameter, relative positron concentration and matching normalized 
propagation-speed range, M, are shown in Fig. 3 for both ISWs and IPWs. As 
it is clearly remarked, only compressive ISWs propagate in this plasma. It is 
observed from Figs. 3(a) and 3(b) that for ISWs/IPWs the increase in the rel- 
ative positron concentration decreases/increases the potential depth and width 
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for other fixed parameters. The same rule applies also to Figs. 3(c) and 3(d), 
which show the variations of potential width/depth (wave amplitude/ width) for 
ISWs and IPWs with respect to different relative matching-speed values, M. 

Figure 4 plots the variations of electrostatic potential (density soliton) profiles 
with respect to the fractional positron concentration, a, and the relativistic de- 
generacy parameter, Rq. From Fig. 4(a) it is clearly observed that, the increase 
in the values of relative positron concentration, a, leads to increase/decrease 
of amplitude/ width of solitary nonlinear structures. On the other hand, from 
Fig. 4(b) it is revealed that, the increase in the values of relativistic degeneracy 
parameter, a, leads to decrease/increases of amplitude/ width of solitary non- 
linear structures. It is therefore confirmed that the narrower ISWs the taller 
it is, which is in accord with the nonlinear theory of soliton in agreement with 
the relation 0qA = est., where, <f) and A denote the amplitude and the width of 
nonlinear waves, respectively. 

Furthermore, Fig. 5 shows the stability regions for nonlinear ISWs and IPWs for two 
extreme cases of mass density. The lower R$ = 1 value corresponds to the mass-density of 
a typical white-dwarf, while, the higher value of i?o = 50 which results in a mass-density 
much above that of a white dwarf, is given to show the effect of higher relativistic degeneracy 
effects on nonlinear wave dynamics, that is, to consider the ultra-relativistic limit of Eq. ([3]) 
(Ro — > oo). It is observed that the stability region for p ~ 1.97 x 10 6 (gr /cm 3 ) lies entirely 
in subsonic, v g < c sr , limit, while for p ~ 2.5 x 10 u (gr/cm 3 ) a large portion of ISWs region 
extends to supper-sonic region, v g > c sr . It should be distinguished that, the speed, c sr , is a 
relativistic speed, despite the name of sound speed and is much higher (c sr ~ 8 x 10 8 (cm/s) 
comparable to the Fermi-speed of electrons in typical solids) compared to that of the ordinary 
plasmas. It is also remarked that, the inclusion of pair modifies greatly the high positron- 
concentration portion of the stability regions (e.g. compare Figs. 2 and 4), as expected. 
Note that, in calculations for the mass density, p ~ 2m p n e o(l — a), there is a a-dependent 
term which accounts for the extra electrons produced/annihilated in the plasma. Since, the 
current analysis is applicable to the wide range of the relativity parameter i.e. the wide 
plasma mass-density range; it can also be applied to low mass-density degenerate plasma 
produced in laboratory. 
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V. CONCLUSIONS 



The Sagdeev pseudo-potential approach was used to investigate the propagation of non- 
linear ion waves in a relativistically degenerate pair plasma in the framework of quantum 
hydrodynamics model. The matching criteria of existence of such solitary excitations in 
terms of the relativistic degeneracy parameter, relative positron density were numerically 
investigated and the allowed matching normalized propagation-speed range, M, for propa- 
gation of such waves was evaluated numerically. It was shown that the electron relativistic 
degeneracy as well as the relative positron concentration has outstanding effects on nonlinear 
wave dynamics of relativistically degenerate plasmas such as that found in white dwarfs and 
the cores of massive planets. Current findings are completely consistent with the previously 
reported similar plasma cases treated with semiclassical assumptions for non-relativistic and 
ultra-relativistic electrons/positrons. This study can be appropriate for applications in in- 
ertial confinement fusion laboratory research, as well as the study of astrophysical dense 
objects. 



12 



[1] A. A. Vedenov, E. P. Velikhov, and R. Z. Sagdeev, Nuclear Fusion, 1 82(1961) [in Russian]. 

[2] M. Salahuddin, H. Saleem and M. Saddiq, Phys. Rev. E, 66 036407(2005). 

[3] A. Esfandyari-Kalejahi, M. Mehdipoor and M. Akbari-Moghanjoughi, Phys. Plasmas, 16 

052309(2009). 
[4] R. S. Tiwari, Phys. Lett., A, 372 3461(2008). 

[5] R. S. Tiwari, Aditya Kaushik, M. K. Mishra, Phys. Lett., A, 365 335(2007). 
[6] A. Mushtaq and H. A. Shah, Phys. Plasmas, 12 072306(2005). 

[7] A. Esfandyari-Kalejahi, M. Akbari-Moghanjoughi and B. Haddadpour-Khiaban, Phys. Plas- 
mas, 16 102302(2009). 
[8] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 052302(2010). 

[9] Prasanta Chatterjee, Kaushik Roy, S. V. Muniandy, and C. S. Wong, Phys. Plasmas 16, 
112106(2009). 

[10] S. I. Popel, S. V. Vladimirov, and P. K. Shukla, Phys. Plasmas, 2 716(1995). 
[11] Y. N. Nejoh, Aust. J. Phys., 50 309(1997). 

[12] S. Mahmood, A. Mushtaq, H. Saleem, New J. Phys., 5 4680(2003). 

[13] M. A. Mahmood, S. Mahmood, A. M. Reza and H. Saleem, Chin. Phys. Lett. 22 632(2005). 

[14] S. Mahmood and N. Akhtar, Eur. J. Phys. D, 49 217(2008). 

[15] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 072101(2010). 

[16] P. K. Shukla, L. Stenflo, and R. Fedele, Phys. Plasmas, 10, 310(2003). 

[17] V. I. Berezhiani, L. N. Tsintsadze, and P. K. Shukla, J. Plasma Phys., 48, 139(1992); Physica 

Scripta, 46 55(1992). 
[18] V. I. Berezhiani, M. Y. El-Ashry, U. A. Mofiz, Phys. Rev., E50 448(1994). 
[19] V. I. Berezhiani, and S. M. Mahajan, Phys. Rev. Lett., 73 1110(1994). 
[20] F. B. Rizzato, Plasma Phys., 40 289(1988). 

[21] H. R. Miller and P. J. Witta, Active Galactic Nuclei (Berlin: Springer) (1987), p. 202 
[22] F. C. Michel, Rev. Mod. Phys. 54 1(1982). 

[23] L. O. Silva, R. Bingham, J.M. Dawson, J. T.Mendonca, and P. K. Shukla, Phys. Rev. Lett. 
83 2703(1999). 

[24] P. Goldreich, W. H. Julian, Astrophys. J., 157, 869(1969). 



13 



[25] M. J. Rees, The Very Early Universe, (Cambridge : Cambridge University Press), (1983). 
[26] C. M. Surko, M. Levelhal, W. S. Crane, A. Passne and F. Wysocki, Rev. Sci. Instrum. 57 
1862(1986). 

[27] C. M. Surko, T. Murphy, Phys. Fluids, B2 1372(1990). 

[28] R. G. Greeves and C. M. Surko, Phys. Rev. Lett. 75 3847(1995). 

[29] V. I. Berezhiani, D. D. Tskhakaya and P. K. Shukla, Phys. Rev. A, 46 6608(1992). 

[30] S. Son and N.J. Fisch, Phys. Lett. A, 329, 76(2004). 

[31] S. Son and N.J. Fisch, Phys. Rev. Lett., 95, 225002(2005). 

[32] M. Akbari-Moghanjoughi, Phys. Plasmas, 17 082315(2010). 

[33] Shapiro, S. L. and Teukolsky, S. A., Black Holes, White Dwarfs and Neutron Stars, (Wiley, 
New York), (1983). 

[34] H. Haug, S. W. Koch, "Quantum theory of the optical and electronic properties of semicon- 
ductors" , World Scientific, 2004, 

[35] Shukla P K, Eliasson B, "Nonlinear aspects of quantum plasma physics" Phys. Usp. 51 
53(2010). 

[36] M. Bonitz, D. Semkat, A. Filinov, V. Golubnychyi, D. Kremp, D.O. Gericke, M.S. Murillo, 

V. Filinov, V. Fortov, W. Hoyer, J. Phys. A, 36 5921(2003). 
[37] S. Chandrasekhar, "An Introduction to the Study of Stellar Structure", Chicago, 111. , (The 

University of Chicago press), (1939). 
[38] L. D. Landau and E. M. Lifshitz, Statistical Physics (Pergamon, London, 1958) 
[39] S. Chandrasekhare, Science, 226 4674(1984). 

[40] W. Misner, K.S. Throne, J.A. Wheeler, Gravitation, Freeman, San Francisco, (1973)763. 

[41] R. Sabry, W. M. Moslem, F. Haas, S. Ali, and P. K. Shukla, Phys. Plasmas, 15, 122308(2008). 

[42] A. P. Misra, N. K. Ghosh and C. Bhowmik, Eur. Phys. J. D, 49 373(2008). 

[43] A. P. Misra, P. K. Shukla, C. Bhowmik, Phys. Plasmas, 14, 082309(2007). 

[44] U. M. Abdelsalam, W. M. Muslem and P. K. Shukla, Phys. Lett. A, 372 4057(2008). 

[45] M. Akbari-Moghanjoughi, Phys. Plasmas, 17, 092304(2010). 

[46] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78 591(2006). 

[47] D. Shaikh and P. K. Shukla, New J. Phys. 10 083007(2008). 

[48] G. Brodin, M. Marklund, and G. Manfredi, Phys. Rev. Lett. 100 175001(2008). 

[49] C. Gardner, SIAM, J. Appl. Math. 54 409(1994). 

14 



[50] A. E. Dubinov, Plasma Phys. Rep., 35, 991(2009). 

[51] A. M. Berthomier, R. Pottelette and R. A. Treumann, Phys. Plasmas, 6 467(1999). 

[52] R. A. Treumann and W. Baumjohann, Advanced Space Plasma Physics , (Imperial Colleage 

Press, London), (1997), page 260. 
[53] M. Akbari-Moghanjoughi, Astrophys. Space Sci., (to be published), DOI10.1007/sl0509-010- 

0499-3. 



15 



FIGURE CAPTIONS 



Figure- 1 

(Color online) Variation of the stability regions of arbitrary-amplitude ISWs/IPWs is 
shown in M-Rq plane for a relativistically degenerate electron-positron-ion plasma with 
respect to different relative positron concentrations, a. The region to the right (yellow) 
correspond to the stability of ISWs and to the left (blue) region to IPWs. 

Figure-2 

(Color online) The variation of Sagdeev pseudo-potential profile with respect to different 
relativistic degeneracy parameters, Ro, relative positron concentrations a and the matching 
normalized propagation-speed range, M, for both ISWs (left column in figure) and IPWs 
(right column in figure). 

Figure-3 

(Color online) Variation of the stability regions of arbitrary-amplitude ISWs/IPWs is 
shown in M-a plane for a relativistically degenerate electron-positron-ion plasma for two 
extreme relativistic cases with respect to relativistic degeneracy parameter, R . The region 
to the right (yellow) correspond to the stability of ISWs and to the left (blue) region to 
IPWs. 

Figure-4 

(Color online) (a) The variation of arbitrary-amplitude ion solitary wave profile is shown 
with respect to different relative positron concentration, a and other fixed plasma param- 
eters, (a) The variation of arbitrary-amplitude ion solitary wave profile is depicted with 
respect to different relativistic degeneracy parameter, R and other fixed plasma parame- 
ters. The increase in thickness of curves in each plot represents the increases in the varied 
plasma parameter. 

Figure-5 

(Color online) Variation of the stability regions of arbitrary-amplitude ISWs/IPWs is 
shown in M-a plane for a relativistically degenerate electron-positron-ion plasma for two 
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extreme white dwarf mass-densities. The parameter p is the corresponding mass-density of 
a typical white dwarf. The region to the right (yellow) correspond to the stability of ISWs 
and to the left (blue) region to IPWs. 
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